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SUMMARY 

This paper presents the structural and aerodynamic aspects of a 
general study of advanced transonic transports. In this study, aircraft 
designed to cruise at Mach numbers of 0.90, 0.98, and 1.15 were comparatively 
analyzed. The wings of all three aircraft employ supercritical sections, 
and the two aircraft with the highest cruise Mach numbers also employ fuselage 
area ruling. 

Structural /aerodynamic characteristics and interactions are investigated 
both parametrically and with the aid of an automated configuration optimization 
program. The effects of replacing conventional aluminum airframe structure 
by advanced filamentary composite (carbon/ epoxy) structure receive particular 
attentioi . The methods employed in the structural /aerodynamic analysis 
are presented in the appendices. 

INTRODUCTION 

It is well known that the commercial transport aircraft industry is 
a vital part of the nation's economy. In order to maintain the present 
U. S. dominance in this field in the face of ril ing foreign competition, the 
next generation of transport designs must incorporate the latest and best 
of advanced technology. For this reason, a study was initiated to examine 
commercial transport aircraft, designed for introduction in approximately 
1980, for the primary purpose of determining the advanced technology areas 
that give the largest improvements in performance and economics over current 
airc .ft. The results of this study are reported in a si'*.iary paper (ref. 1) 
and papers in the specialty areas of aerodynamics and structures, propulsion 
(ref. 2), and economics (ref. 3). This paper reports the methods used 
and the results obtained in the areas of aerodynamics and structures. 
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The study is limited to transcontinental aircraft which would not 
produce a sonic boom on the ground. The aircraft studied are conventional 
swept wing/fuselage configurations designed for a range of 2700 n. ml. with 
a payload of 200 passengers and no cargo. 

One effect investigated was that of Increasing cruise speed. Increased 
cruise speed is valuable, not only because of a possible competitive market 
advantage, but also because of the resulting Increased aircraft productivity. 
Aircraft were studied that are designed to cruise at Mach numbers of 0.90, 

0.98, and 1.15. For the purposes of this study the aircraft are denoted 
respectively as CVT (conventional transport), ATT (advanced technology transport), 
and TST (transonic transport). Cruise at Mach number 0.90 represents approximately 
the highest cruise speed achievable without significant wave drag for an 
aircraft using the supercritical wing and no fuselage area ruling. Cruise 
at Mach number 0.98 represents approximately the highest cruise speed achievable 
without wave drag by using the supercritical wing and fuselage area ruling. 

A cruise speed of Mach number 1.15 is of Interest because it represents 
approximately an upper limit for cruise with no sonic boom reaching the 
ground. Because temperatures at cruise altitude are lower than those at sea 
level, the speed of sound is lower and as a result an aircraft flying at 
Mach 1.15 at 40,000 feet Is actually flying at a velocity that corresponds to 
Mach 1 at sea level. Both theory and flight tests Indicate that no boom 
is heard on the ground under these conditions. Recent results suggest that 
Mach 1.15 may be too high; Mach 1.08 may be more reasonable to allow for terrain 
and weather fluctuations encountered In flying across the country. A cruise 
Mach number of 1.08 at 40,000 feet would correspond to a sea level Mach number 
of 0.94. For the purposes of this study, it is felt that the results pertaining 
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to the aircraft designed to cruise at Mach 1.15 would be very close to 
those for an aircraft designed to cruise at Mach 1.08. 

The aircraft of this study were analyzed with the aid of TRANSYN-TST, 
an aircraft synthesis program developed by ACMD personnel. The TRANSYN-TST 
computer program consists of a control program and a number of subroutines 
to do the various specialty tasks required to synthesize an aircraft design. 
Comparisons were facilitated by using the same computing methods and ground 
rules (range and payload requirements) for all aircraft. The methods used 
in the aerodynamic and structures subroutines of TRANSYN-TST are described 
in detail in the appendices. 

The studies of advanced aerodynamics technology involved application 
of the supercritical wing, which was developed by Dr. R. Whitcomb at Langley 
Research Center. With this wing it is possible to delay the onset of wave 
drag and/or to use thicker wing sections with lower wing sweep. Proper 
fuselage area ruling was used to minimize the wave drag for the aircraft 
designed to cruise at Mach 1.15 and to eliminate it for the aircraft designed 
to cruise at Mach 0.98. In the structures/materials area, the effect of 
advanced composite materials on aircraft structural weight and performance 
was investigated. In particular, the effect on gross takeoff weight of 
using carbon/epoxy composite material in the primary load carrying structure 
was analyzed. Combining advanced structures technology with advanced 
aerodynamics leads to interesting tradeoffs in aircraft design. For example, 
more efficient structures allow the aircraft to be designed to achieve 
higher aerodynamic efficiency. These structures and aerodynamic interactions 
are described in this report. 



- 4 - 


METHODS OF ANALYSIS 
Structures 

Traditionally, preliminary weight estimates of aircraft are made 
using empirical methods based on the weights of existing aircraft. These 
methods, however, are undesirable for studies of the type reported here 
for two reasons. First, since the weight estimating formulas are based 
on existing conventional aircraft, their application to unconventional 
configurations (e.g., those with area ruled bodies) is suspect. Secondly, 
they provide no straightforward manner for assessing the impact of advanced 
technology structures and materials (e.g., filamentary composite materials). 
For the present study the weights of the load carrying wing and fuselage 
structures were determined by analytically based methods. Since it was 
necessary to be able to rapidly evaluate a large number of specific designs, 
the methods employed were based on idealized vehicle models and simplified 
structural analysis. 

The configurations considered in this study are of the conventional 
wing-body type. As shown schematically in figure 1, the fuselage structure 
is an integrally stiffened shell stabilized by frames and the structure is 
sized by longitudinal bending. The wing is a mur;i-web box beam designed 
by span-wise bending and shear. The lower portion of the wing-fuselage 
carry-through structure carries the span-wise bending, shear, and torsion 
loads introduced by the outboard portion of the wing; the entire carry- 
through structure carries longitudinal bending loads. 

The method of determining the fuselage structural weight is based on 
simple beam theory and is described in Appendix A. The bending moment 
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distribution resulting from a 2.5 g maneuver is computed first. The 
analysis then proceeds on a station-wise basis, calculating the minimum 
amount of material at each longitudinal station to satisfy static strength 
and stability requirements. The fuselage shell is sized by buckling with 
a compression strength cut-off. In addition, the shell is checked for 
minimum gage and tensile strength requirements and resized if necessary. 

For the buckling analysis, the integrally stiffened shell is analyzed 
as a wide column and the frames are sized by the Shanley criterion. 

The wing weight analysis is described in Appendix B. The weight 
of the structural wing box is determined by calculating the minimum amount 
of material required to satisfy static buckling requirements at a series 
of span-wise stations. The covers of the multi -web box are sized by 
buckling due to local instability and the webs by flexure- induced crushing. 
Required shear material is computed independently of buckling material. 
Aeroelastic effects are not accounted for directly. For the carry-through 
structure, buckling, shear, and torsion material are computed independently 
and summed. 

The weight analyses of the fuselage and wing described above do not 
completely account for the total weight of these structures. Specifically, 
weight penalties due to fasteners, cutouts, surface attachments, uniform 
gage requirements, manufacturing constraints, aeroelastic effects, etc., have 
not been directly determined. These weight items, conmonly aggregated in 
a term called "non-optimum" weight, are determined by an empirical method 
based on existing aircraft, as described in Appendix A. The resulting 
"non-optimum factors," are 2.9 for the fuselage and 2.6 for the wing. 

The weights of all the other items making up the dry weight of 
the stuc(y aircraft are estimated by empirical formulas derived from 
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regression analysis of existing aircraft. Such formulas predict the 
weight of these items satisfactorily for comparative studies. 

The structural and weight analysis methods were used to estimate the 
weights of two existing aircraft, the 707-120B and the 747B, and the results 
are compared with the actual weights in figure 2. These aircraft were 
not used in the regression analyses to determine the non-optimum factors 
and therefore the comparison provides a good check on the methods. The 
agreement is considered excellent for the 707-120B and good for the 747B. 

The descrepancy in the case of the 747B is primarily in the fuselage weight, 
and may be because the 747B airframe is a more recent design then the airframes 
of the 707-1 20B and the aircraft used in the regression analyses. 

Aerodynami cs 

The methods used to estimate the lift and drag characteristics of 
the study aircraft are described in detail in Appendix D and sumnarized 
in figure 3. The major differences between the study aircraft configurations 
are shown in figure 4. These factors determined the aerodynamic methods 
which were applied. The wing sweep angles chosen for each configuration 
resulted from constraining the Mach number perpendicular to the wing quarter 
chord line to 0.74. This constraint and a constraint of 10% on the wing 
thickness ratio normal to the wing quarter chord were applied to insure 
good wing section performance. These aerodynamic constraints are shown 
as a function of cruise Mach number in figure 5. 

In order to compare the relative landing qualities of the study aircraft, 
the performance of each aircraft with hirh-lift devices was estimated. 

A computer program based on the relationships described in reference 4 was 
used to calculate the effectiveness of high-lift devices on various wings 
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with different wing sweep, aspect ratio, taper ratio, and thickness ratio. 

All configurations were analyzed with the same type of high-lift devices, 
double slotted trailing edge Fowler flaps and leading edge slats. The 
aircraft approach speed was computed at 1.3 times the aircraft s '11 speed, 
based on the estimated maximum lift coeffic ^nt ( c i_max^ dnc * an aircraft 
maximum landing weight with 0.7 of the mission fuel consumed. 

CHARACTERISTICS OF NOMINAL VEHICLES 

The nominal aircraft configurations are shown in figure 6. All the 

aircraft were designed for a range of 2700 nautical miles with a pcyload 

of 200 passengers and no cargo. The nominal values of aspect ratio, wing 

2 

loading, and engine bypass ratio are 7,120 lbs/ft , and 2, respectively. 

The aircraft are drawn to scale to show the relative sizes and the effect 
of fuselage area ruling. The gross takeoff and operating empty weights 
are 235,000 lbs and 130,000 lbs for the CVT, 250,000 lbs and 140,000 lbs 
for the ATT, and 600,000 lbs and 235,000 lbs for the TST. Thus, due primarily 
to wave drag, the TST is significantly heavier than the other two con- 
figurations. 

The principal difference in the fuselage structure of the configurations 
considered here ana those of current subsonic jet transports is the presence 
of fuselage area ruling. As mentioned earlier, the M 1.15 vehicle (TST) 
and the M 0.98 vehicle (ATT) are area ruled, whereas the M 0.90 vehicle 
(CVT) is not. Since the TST is the most highly area ruled of the configura- 
tions, this configuration will be used to illustrate the effects of area 
ruling on body structural weight. 

The wave drag encountered by the TST results in a large fuel requirement 


which leads to a high gross takeoff weight (nearly three times that of the 
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other study configurations). Aerodynamic support of this high weight requires 
a large wing. Thus the fuselage of this configuration is relatively highly 
loaded and severely area ruled when compared to *he CVT and the ATT. Figure 
7a shows the body contour; the area ruling is clearly seen. Figure 7b 
shows the distribution of longitudinal bending moment on the body. Note 
that this distribution has a sharp peak which occurs over the area ruled 
portion of the body. Since body shell weight per ur.it. surface area* which 
is shown in figure 7c, increases with increasing bending moment and decreases 
with increasing radius, this weight is large in the area ruled portion 
of the body, especially aft of the wing carry-through structure. (Recall 
that this structure helps to carry the longitudinal bending loads.) The 
bump in the unit weight at the aft end of the body is due to introduction 
of the tail load. However, the effect of the increase in unit weight at 
this body station on the overall fuselage structure weight is small because 
of the relatively small fuselage diameters in this region. The shell is 
sized by the minimum gage constraint at the front of the fuselage and by 
buckling elsewhere, except ~ : n the region of the peak bending moments where 
the compression strength cutoff is active. 

The principal differences in the wing structure of the configurations 
considered here and those of current subsonic transports are higher sweep, 
higher* aspect ratio, and thicker supercritical wing airfoil sections. 

The higher sweep and aspect ratio result in relatively higher wing weights 
for the study configurations in spite of the greater wing thickness. The 
“optimized" configurations employing advanced composite materials have 
particularly high aspect ratios, significantly higher than existing 
transport aircraft. Since aeroelastic effects have not been accounted 
for directly (they have been accounted for Indirectly for moderate aspect 
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ratios through the "non-optimum" factor), the wing weight estimates for 
these configurations must be regarded as possibly optimistic. 

SENSITIVITIES OF NOMINAL VEHICLES TO CONFIGURATION PARAMETERS 

The sensitive ies of the previously discussed nominal configurations 
will be presented in the form of variations in airframe weights s (L/D)^, 
and gross takeoff weight (Wq T q) with changes in the configuration parameters 
aspect ratio (M), wing loading (W/S), thickness ratio (t/c), and sweep 
(a). In each variation, the parameters not being changed are held constant at 
their nominal values. The payload-range characteristics of the vehicles 
are also held constant at their nominal values (200 passengers, 2700 n. mi.). 
Operating weight empty fraction (0 WE/W GT q) a*.- (L./D)^ a^e used &s measures 
of aircraft structural efficiency and aerodynamic efficiency, respectively. 

It should be noted, however, that because of the resizing process a large 
change in OWE will result in only a small change in OWE/W^. Configuration 
gross weight, which is closely related to direct operating cost, is used 
to indicate overall aircraft efficiency. 

The sc 'sltlvities of the wing, fuselage, and OWE weight fractions 
to A*s are shown in figure 8a. The win^ weight fraction increases with 
increasing & due to increasing span; the fuselage weight fraction increases 
due to decreasing root chord and also, for the ATT and TST, due to increas- 
ingly severe fuselage area ruling. The sharp uptrend for the TST above 
& a 7 is a size effect in which attempts to accomplish the mission by 
increasing the aircraft gross weight are uverualanced by increasing structural 
weight. Above & * 8, *he TST cannot carry 200 pa"*engers 2700 n. mi. at 
any gross weight. The effect for all configurations is increased (WE 



fraction with Increasing A^. However, also has a strong effect on 
(L/0 )max as shown In figure 8b, and the increasing aerodynamic efficiency 
with increasing M (due to reduced induced drag) tends to balance the de- 
creasing structural efficiency, with the result that W~ TO is relatively 
insensitive to & for the CVT and the ATT. The TST, being marginally able 
to perform the mission, is sensitive to almost any parameter, and the TST 
gross weight exhibits a sharp minimum at about the nominal value. 

The wing weight fraction, figure 9a, decreases with increasing W/S 
due to decreasing wing area. The fuselage weight fraction varies only 
slightly due to changing root chord length and fuselage area ruling, except 
for the TST. The net effect is increasing structural efficiency (decreasing 
OWE weight fraction) with increasing W/S. However, since (L/D)^ decreases 
with increasing W/S, the situation is the same as for & ; the structural 
and aerodynamic effects tend to neutralize each other and Uqjq is relatively 
insensitive to l'/S (figure 9b). The exception is the TST which improves 
with higher W/S. This trend, however, is constrained by landing speed 
limitations; the effect on landing quality is considered later in this 
paper when “optimized" configurations are discussed . 

Figure 10a shows that wing weight fraction increases slightly with 
decreased t/c. The fuselage weight fraction decreases slightly for the 
area ruled ATT and TST configurations because a thinner wing results in 
a less severely area-ruled body. The OWE fraction is relatively insensitive 
to t/c. Since ( *-/D)^ x is also very insensitive to t/c withir reasonable 
limits, as shown in figure l^u, this parameter is relatively unimportant 
for configuration optimizatio n a r - ’ was held at the aerodynamically 

constrained upper limit for ea ** ..drfon throughout the study. 

Figure 11a sh ows that wing weight fraction Increases with Increasing 
A while fuselage weight fraction remains fairly constant. The TST again 
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is much more sensitive than the other configurations. Since the lowest 
Wqtq occurred at the aerodynamically constrained minimum sweep value, 
as shewn in figure lib, the sweep was held at this minimum value for each 
configuration throughout the study. 

ADVANCED COMPOSITE MATERIALS 

Because tne best prospect for substantially improving the performance 
of aircraft appears to be the use of advanced composite materials, the 
potential impact on transport aircraft performance of these materials 
was investigated in this study. Two advanced composites materials were 
considered: boron filaments in an epoxy matrix (B/E) and carbon filaments 

in an epoxy matrix (C/E). For the purposes of this study, it was assumed 
that these materials are layed up in isotropic laminates consisting of 
orthotropic lamina, each lamina being composed of unidirectional filaments 
in a homogeneous matrix. The manner in which the material properties of 
the composites were characterized is derived and discussed in Appendix C. 
Only complete substitution of composite materials in both the wing and 
fuselage load carrying structures was considered in the present study. 

It should be noted that selective substitution may be more cost effective 
than total substitution and that composite materials may also be attractive 
for use in other structures such as tails and landing gear. 

Before investigating the effect of composite materials on the study 
configurations, it is instructive tc consider the material properties as 
given in Table Cl of Appendix C. In the last column, p/*1T is tabulated 
with respect to aluminum alloy; this quantity is a good Indicator of the 


P 
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weight of buckling critical structures. The table indicates that a composite, 
buckling critical structure will be about 60% of the weight of an equivalent 
aluminum alloy structure. Since the weight of strength critical structure 
is proportional to p/F, the table indicates that strength critical structure 
will show an even greater weight reduction with the use of advanced composites 
than will buckling critical structure. 

The weights which result when the wing and fuselage structures are 
made entirely of C/E are shown in Table 1 and compared with the weights 
of the same size aluminum airframes. Since these configurations were 
analyzed at the same gross takeoff weight, the aircraft which results in 
lower operating weight empty (OWE) will be able to carry more fuel and 
thus have longer range with the same payload. The net reductions in 
wing and fuselage weights for the three configurations range from 33% 
to 43%. This translates into reductions of from 15-1/2% to 19-1/2% 
in the OWE. The reduction is largest for the TST because this vehicle 
is heavier and more highly loaded than the other two vehicles and thus 
is more strength critical. 

Rather than hold gross takeoff weight constant during the material 
substitution, it is more interesting to see the effect of the substitution 
when holding the payload and range fixed. In this substitution, the 
magnification effect (lower structural weight means lower fuel weight means 
lower structural weight, etc.; is accounted for. Figure 12 shows that 
composites reduce the gross weight of the nominal CVT slightly (15%), the 
nominal ATT slightly more (19%) and the nominal TST substantially (48%). 

In every case, the C/E composite is better (lighter weight) than the B/E, 
as is also obvious from Table Cl. Therefore, further analysis of C/E was 
emphasized in the present study. The gross takeoff weights of the CVT 

# 

I 



and ATT are signficantly les* than that of the TST, even when the latter 
is made of composite material. Since the airframe cost of an aircraft 
made from advanced comoos ite materials will be different than the cost 
of a conventional aluminu. alloy aircraft, material selection must be 
based on overall cost effectiveness, and not solely on the basis of least 
weight. It is clear, however, that the TST is a much more attractive vehicle 
when made of C/E than when made of aluminum. 

The sensitivities of the gross takeoff weights (for the specified 
mission of 2700 n. mi. with 200 passengers) of the C/E configurations 
to the configuration parameters ^ , W/S, t/c, A are shown on figures 13-16. 
Also shown for reference are the sensitivities of the aluminum configurations 
discussed earlier. Sensitivity of gross takeoff weight to M , figure 13, 
shows the same affect for C/E as for aluminum in the case of the CVT and 
ATT. For the TST, the C/E configuration is much less sensitive to PK. 

The (L/D)^ is less for the lighter C/E configurations than for the corre- 
sponding aluminum configurations due to the smaller relative wing area 
and fixed fuselage size (wing loading is held constant). Substitution 
of C/E for aluminum results in little change in sensitivity to W/S, figure 
14. The C/E TST, however, is much less sensitive to W/S than is the aluminum 
TST, but least weight still results at higher W/S than nominal. The material 
substitution has the same affect oi the sensitivities to t/c, figure 15, 
and to A, figure 16, as it had on the sensitivities to and W/S; the CVT 
and ATT remain relatively insensitive while the C/E TST becomes much more 
insensitive than the aluminum TST. 
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OPTIMIZED CONFIGURATIONS 

As has been shown, the vehicles of this study have strong structural/ 
aerodynamic interactions and in many cases structural efficiency may be 
traded for aerodynamic efficiency. To investigate these interactions, 
the synthesis program was coupled with a parameter optimization program 
called AESOP (refs. 5 and 6). AESOP adjusted the values of A, W/S, and 
engine bypass-ratio to obtain minimum Wqjq. The parametric results presented 
earlier show that t/c and A are best at their constrained values and therefore 
these parameters were not allowed to vary in the optimization. Gross takeoff 
weight was chosen as the quantity to be minimized since it also closely 
reflects the direct operating cost (unless materials are changed). Since 
some unconstrained configurations optimized at high U/S and low & which gave 
high landing speeds, the affect of a performance penalty for landing speed 
was also investigated. Consideration of optimized configurations is particularly 
important when evaluating aircraft made of composite structure because 
the optimum configurations of these aircraft may be considerably different 
than for aircraft made of aluminum. Thus, the resultant weight advantage 
of composite materials might be considerably greater than that predicted 
by substitution without configuration reoptimization. 

Six configurations were optimized in the manner just described. 

These configurations utilized the two materials aluminum (Al) and C/E 
for the three study vehicles CVT, ATT, and TST. The weight breakdowns 
of the configurations are shown on figure 17; also shown for conparison 
are the weights of the nominal configurations. The gross takeoff weights 
of the CVT and ATT, Al and C/E configurations have been reduced by about 
5 % relative to the nominal configurations. The weight of the Al TST has 
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been reduced by about 15%. However, the weight of the C/E TST has not 
been reduced relative to the nominal due to the fact that the wing size 
of this configuration is limited by a minimum fuel volume constraint. It 
is significant that the relative ranking in terms of gross weight of the 
configurations is the same for the optimized configurations as for the 
nominals. As was expected, gross weight is insensitive to the configuration 
parameters in the vicinity of the optimums. 

Drawings and some of the important parameters of the three optimized 
A1 and C/E aircraft are shown in figures 18 and 19, respectively. In redesigning 
the aircraft with C/E material, the increased structural efficiency allowed 
increased aerodynamic efficiency. Optimization of the aircraft configuration 
for each case insured that the best balance between weight and aerodynamics 
was obtained. The largest reduction in aircraft gross takeoff weight due 
to the use of C/E material was for the TST configuration. As mentioned 
previously, this is partly because the A1 TST is operating close to its 
ultimate mission capability. 

For all configurations the use of C/E results in an increase in the 
optimum wing aspect ratio. The aspect ratio of the CVT increased from 
6.8 to 9.1 with C/E, that of the ATT increased from 8.1 to 12.4, and that 
of the TST went from 6.4 to an aspect ratio of 9.8. At the same time, 
the optimum W/S decreased in all cases: from 123 to 113 psf for the CVT, 

123 to 106 psf for the ATT, and 127 co 100 psf for the C/E TST compared 
with the A1 TST. Since the increased structural efficiency results In 
large M for the C/E configurations, questions as to possible flutter 
problems arise. As stated previously, the present TRANSYN-TST synthesis 
program does not contain any flutter calculation, and no additional structural 
weight is added to increase the wing stiffness to avoid flutter. 
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Because of the lower wing loadings and higher aspect ratios of the 
optimized C/E configurations, their approach speeds are lower than for 
the optimized A1 configurations. The effect Is strongest in the case of 
the TST, where the approach speed has decreased from 190 knots for the 
A1 TST to 168 knots for the C/E TST. The effects of W/S, %, and A on 
the C|_|^ x and the resultant effect on approach speed for the three nominal 
aluminum configurations are shown in figures 20, 21 and 22. As mentioned 
previously, the approach speeds of all configurations were estimated for 
double slotted flaps and leading edge slats. 

Using the parameter optimizer AESOP and including a performance penalty 
for high approach speed results in a change in the optimized A1 ATT as 
shown in figure 23. On the left of the figure is the ATT configuration 
optimized without a penalty for high landing speed and on the right is 
the ATT configuration optimized with the landing speed penalty. The result 
of including this penalty is to increase the A and decrease the W/S. These 
design changes result in a slightly higher Wqjq because of the higher operating 
weight empty due to the increased wing and fuselage structural weight. 

The decrease in approach speed is substantial, from 170 knots for the uncon- 
strained A1 ATT to 151 knots for the constrained A1 ATT. Figure 24 shows 
the affect of an approach speed constraint on the optimized C/E TST con- 
figuration. As with the ATT, the inclusion of this penalty for high landing 
speed results in higher A and lower W/S and leads to higher operating 
weight empty and higher Wqjq. The decrease in approach speed is less for 
the TST than the ATT, having dropped from 168 knots to 158 knots. Because 
of the large penalty in overall performance of this small landing speed 
reduction, possible consideration should be given to other means of obtaining 
high-lift at takeoff and landing, such as variable geometry. 
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The structural /aerodynamic characteristics of advanced transonic 
transports have been investigated. Specifically, aircraft designed to 
cruise at Mach numbers of 0.90 (CVT), 0.98 (ATT), and 1.15 (TSS) were 
studied. It was found that the CVT and ATT designs had similar ch"-acter- 
istics and performance (as measured by gross takeoff weight) for the same 
payload/range requirements (200 passengers, 2700 n. ml.). The principal 
difference between the CVT and the ATT is the use of area ruling in the 
latter, resulting in a slight increase in fuselage structural weight. The 
gross takeoff weight of the TST, however, was found to be more than twice 
that of the other two designs for the same mission if all designs employ 
conventional aluminum airframes. The reason for the poor performance of 
the TST is the large wave drag of this configuration. 

Parametric analyses were made for the three designs. These analyses 
show that both aspect ratio and wing loading had strong effects on aero- 
dynamic efficiency as measured by (L/D)^ and structural efficiency as 
measured by operating weight empty. However, these effects tended to 
cancel each other and the net effect of changes in either parameter on 
the gross takeoff weight was small. The parametric results also showed 
that the values of wing sweep and thlckness-to-chord ratio should be held 
at their aerodynamically constrained values. 

Most important of the advanced technologies investigated was that of 
advanced composite materials. Substitution of carbon/epo*y material for 
aluminum in both the wing and fuselage primary structure (without resizing 
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the configurations) resulted in 35-40% reductions in the weights of these 
structures. This translated into 15-20% reductions in operating weight 
empty. If the composite designs are resized to the same payload/range 
criteria as the aluminum ones, the change to composite material results 
in decreases in gross takeoff weight of 15% for the CVT, 19% for the ATT, 
and 48% for the TST. 

Six configurations (consisting of the three designs CVT, ATT, and 
TST made of the two materials aluminum and carbon/epo*y) were "optimized" 
with the aid of a parameter optimization program. In this optimization, 
the values of aspect ratio, wing loading, and engine bypass -ratio were 
altered to obtain minimum gross takeoff weight, both with and without 
an approach speed constraint. As expected, the gross takeoff weight is 
insensitive to variations in the configuration parameters in the vicinity 
of the optimums. For all three vehicle designs the use of carbon/epoxy 
as opposed to aluminum results in a higher value of the optimum wing aspect 
ratio due to increased structural efficiency being traded for increased 
aerodynamic efficiency. Introduction of an approach speed constraint results 
in modest performance penalties for the ’ T configurations but, due to high 
wing sweep, severe penalties for the TST configurations. 
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APPENDIX A - BODY WEIGHT ANALYSIS 

In this appendix, the methods used to estimate the body weights 
of the study configurations are presented. These methods are the 
same as those presented in reference 7, and will only be given 
briefly here. Only configurations with circular cross-sections are 
considered. 

The first sttp in determining body wetght is computation of the 
longitudinal bending moment distribution. To this end the longitudinal 
weight and lift distributions and vehicle center of gravity are 
computed using the body profile, R(X), determined as described f • 
Appendix E, and the wing (see Appendix B), tail, and propulsion j'tem 
weights. The longitudinal bending moment distribution is then obtained 
by considering vehicle loading due to a static 2.5g pull up maneuver 
(factor of safety of l.b). See reference 8. 

The bending moment M is used to compute at each longitudinal body 
station, x, the axial stress resultants at the point of maximum stress 
due to bending. These are, using the flexure formula for beams,* 




+ N 


a 

x 


N 


+ 

x 





a 


for compression and for tension, respectively, where I' is the moment 
of inertia of the body snell divided by the thickness and N x a is the force 
per inch due to axial acceleration. Note that pressure loads are not 
used to relieve compressive stresses. The value of gage pressure, 

P g = 8.77 psi, was held throughout the study. The stress resultant in the 


♦symbols are the same as in reference "i 
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circumferential direction Is 
N y * RP 9 K p 

wiiere Kp accounts for the fact that not all of the shell material is 
avail able for resisting hoop stress. It is assumed that the shell 
material is uniformly distributed around the circumference. 

The equivalent isotropic thicknesses of the shell are given by 





MAX (N*. N y ) 



for designs limited by compression, tension, and minimum gage respec- 
tively where F C y is the compressive strength and is the tensil 
strength. The material is assumed to have elasto-plastic behavior. A 
fourth thickness which must be considered is that for buckling critical 
structures, t«. \ expressions for this quantity will now be derived. 

The nominal vehicles of this study have integrally stiffened shells 
stabilized by ring frames. In the buckling analysis of these structures, 
the shell is analyzed as a wide column and the frames are sized by the 
Shanley criteria (ref. 9). The buckling equation for wide columns 
(ref. 10) gives t c as 
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where d is the frame spacing and E is the modulus of elasticity. Values 
of the shell efficiency e w are given in Table Bl. The Hanley criterion is 
based on the premise that the frames act as elastic supports for the wide 
column; this criterion gives the smeared equivalent isotropic thickness 
of the frames as 


tp = 2 r 2 ~y 

D 


Z~\l _ C F N x 
***** 


where Cp is Shanley : s constant and Kp-j is a frame geometry parameter. See 
reference 7 for a discussion of the applicability of this criterion and for 
a detailed derivation of the equations presented here. If the structure 
is buckling critical, the total thickness is 



Minimizing t with respect to d results in 
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so that the shell is three times as heavy as the frames. 

At aach fuselage station x, the shell may be sized by compression, 
tension, minimum gage, or buckling and the frames are sized by general 
instability. For minimum weight, *t is desired to find the least total 
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thickness t = = tp as a function of d which satisfies the four conditions 

resulting from the combinations of shell and frame criteria. Three of 
these conditions are monotonically decreasing with respect to d and the 
fourth, the buckling condition, has a minimum. The minimum t is obtafned 
by a simple search procedure. The "ideal" or analytical weight of the 
body structure is then obtained by simulation as 

“i “ 2 * P £y ^ AX < 
length 

where quantities subscripted 1 depend on x. 

The "non-optimum" body weight (due to fasteners, cutouts, surface 
attachments, uniform gage penalties, manufacturing constraints, aeroelastlc 
effects, etc.) has yet to be determined. This weight is computed by an 
empirical method. The analysis described above was applied to five 
existing aircraft and the factor was determined which, when multiplied 
by the “Ideal" weight, gave the best fit with the actual weights. This 
factor, called the "non-optimun" factor, was found to be 2.9; thus the 
non-optimun portion of the weight Is nearly twice the analytically 
determined weight. 
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APPENDIX B - WING WEIGHT ANALYSIS 

The lifting planforms of the configurations of the present study 
are tapered, swept wings with straight leading and trailing edges. 

The wing structure was assumed to be a rectangular multi -web box beam 
with the webs running in the direction of the structural span. 

The first step in computing the wing weight is the determination 
of the geometry of the structural wing box. In terms of the input 
parameters Wtq, (W/S) , M , R-^p and the dependent parameters 
wing area, span, root chord, tip chord, and trailing edge wing sweep 
are computed from 

~ W T0 

S " WsT 

b = /RTT 

C R = b (1+R tap ) 

C T = ^TAP C R 

2 Cp 

tan a te = tan a le + — y ( R TAP D 

(See figure Bl). It is assumed that the forward 15% and the aft 15% of 
the streamwise chord are required for controls and high lift devices, 
leaving 70% for the structural box. The intersection of this box with 
the body contours determines the location of the rectangular carry- 
through structure. 

The dimensions of the structural box and of the carrythrough 
structure are now determined. (See figure B2). The structural semi- 
span, b $1 is assumed to be at the quarter-chord line, whose sweep is 
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given by 

3 1 

tan A $ = 4 tan a le + ^ tan Aj E 

Thus, 

b - W r 

k - w 

s “ 2 cos Ag 

The streamwise chord at any point on the wing is given by 

r ^ = C R “ b/2 ^ C R 
In particular, at the wing-body intersection, 

C R C R “ F" ^ C R _ C T^ 

The structural root and tip chords are 

C SR = J C R 
C ST = * 7 C T 

respectively. In terms of y, the structural and total chords are given by 

r sM = C SR " b^ ^ C SR " C ST^ 

r ’ C R " ^ C R " C T^ 

and the thickness of the wing at any spanwise station y is 
t (y) • R t r 

It Is assumed that the fuel Is carried within the wing structure; the 


volume of this structure Is found as follows: 
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V W = 2 J g % r s (y) tM dy = W C SR W C 

- 1 ' 4R t / rbs [ C R - 6^ (C R - V] 2 «y + °- 7 ¥r w C 
"0 

V W - 1 - 4R t b s[ C R- C R< C R- C T) + T( C R- C l) 2 ] +0 - 7R t C R 2w C 

This equation is based on flat upper and lower surfaces and neglects the 
volume of the structure. 

The weight analysis proceeds in a station-wise manner along the 
y (structural semi-span) coordinate. At each station, the lift load, 
center of pressure, inertia load, center of gravity, shear force and 
bending moment are computed. For the inertial load, it is assumed that 

the fuel weight W CT is distributed uniformly with respect to wing volume 

h 1 u 

"FT 

so that the inertial load at y is t f 1 - V(y) where V(y) is the volume 

V W 

outboard of y; this volume has centroid Cg (y) with respect to station y. 
An estimate of the wing structural weight is included in Wpy for this 
calculation but the calculation is not redone when the actual structural 
weight has been computed. It is assumed that the lift load is distrib- 
uted uniformly over the wing planform area so that the lift load at y is 
(W/S) A(y) where A(y) is the area outboard of y; the centroid of this area 
is denoted C p (y). 

The shear force is 

n e 

W 

F.(y) = n K, C(W/S) A — P- V - E h(y.-y) M H 

S S v H 1=1 i 

where n. is the number of engines mounted on the semi span, W. is the 
e 6^ 

XL xL • 

weight of the 1— engine, y^ is the location of the 1— engine, and 
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l.Yi >y 
Myj-y) = l > 

|o, y. < y 

The bending moment Is 

W n e 

M(y) = n Kj [(H/S) A C p - V C g - f h(y r y) W e . (y r y)J 

for a static pitch-plane maneuver of load factor n. 

The dimensions and weight of the structural box can now be calculated. 
Reference 10 Indicates that the critical Instability mode for multiweb 
box beams Is simultaneous buckling of the covers due to local instability 
and of the webs due to flexure induced crushing. This reference gives the 
solidity of least weight multiweb box beams as 



where e and e depend on the cover and web geometries (see Table B1). The 
solidity is defined as 



where Wg END is the weight of bending material per unit span and p is the 
material density. 

The weight of shear material is 

W SHEAR M = ^ 

The optimum web spacing (see figure B2) is computed from 
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The equivalent Isotropic thicknesses of the covers and webs are 


1 



respectively, and the gage thicknesses are 



t = K t u 
9 W 9 W W 


Values of e, e, e r , e r , e u , K , K are found in Table B1 for various 

^ u w "c 

structural concepts (ref. 10). If the wing structural semi-span is 
divided into N equal segments, the total "ideal" weight of the wing box 
structure is 

2 b A 

W B0X = U ^ ^ W BEND 1 + W SHEAR 1 ^ 

The wing carrythrough structure consists of torsion material in 
addition to bending and shear material. The bending material is computed 
in the same manner as that of the box except that only the longitudinal 
component of bending moment contributes, i.e. (let t = t (o) and 
M q = M (o)), 



The weight of the bending material ii then 
W BEND C s p Z C Sr t o W C 
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The quantities d w » t w and t c are computed in the same manner as for the 

box. The weight of the shear material Is (F = F (o) ) 

^o ^ 

W SHEAR C s p of" W C 

The torque on the carryth rough structure is 
T = M q sin A $ 

and the weight of the torsion material is then 

H . pT( V C S R )W C 

w T0RSI0N r 7—f 

t t V G c 

0 SR S 


Finally, the "ideal" weight of the carrythrough structure is computed 
from 

W C = W BEND q + W SHEAR C + W T0RSI0N c 

As is the case with the fuselage structural weight, "non-optimum" 
weight must be added to the ideal weight to obtain the wing structural 
weight. This weight was estimated by applying the above analysis to 
five existing aircraft and determining the factor which, when multiplied 
by the ideal weight, gave the best fit with the actual weights. The 
resulting expression is 

W WINti s 2,6 ( W B0X + W C> 



- 31 - 

APPENDIX C - COMPOSITE MATERIAL CHARACTERIZATION 


The composite materials considered in this paper consist of orthotropic 
lamina formed into quasi -Isotropic (two-dimensional ly isotropic) laminates. 
Each of the lamina is composed of either boron or carbon filaments placed 
uni directionally in an epoxy matrix. Such a quasi-isotropic laminate has 
been shown to result In minimum weight (ref. 11) for the structures of 
the type considered in the present study. Because the lair, ate Is 
two-dimensional ly isotropic, it may be analyzed in the same manner as a 
metal ic material. It remains only +o find the equivalent composite 
properties to be ised in the analysis. 

Because the filament orientations must be dispersed in several 
directions to obtain isotropicity, the uniaxial stiffness of the composite 
will be reduced from that of the individual lamina. This reduction factor 
can be determined using the methods of reference 12. Let a quasi -isotropic 
laminate of total thickness T be In a state of plane stress in its plane 
of isotropicity. Then from (3-36)* the midolane stress resultant-strain 
relations are 


N x * A 11 e x + A 12 V 
N y = A 12 e x + A 22 e y 


where the xy plane is the plane of isotropicity. By analogy with Isotropic 
materials we define the equivalent isotropic modulus of elasticity as 


c = i A,, = =• A 


. 1 
'll " T 


22 


Let the laminate have n layers of equal thickness t * Then from 


(3-? c ' 


A 11 * 


^ ^11 ^k ^ h k~ h k-l ^ = * ®11 ^k 


*all equations refer to Ref. 12. 
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From (2-35), for each layer. 


(Qll) k ~ cos 0 k + 2 (Q^ 2 '^ 2 Qgg) si n O^cos Oj^+Qg^ sin 0^ 

where fi k Is the orientation of the filaments in the kth lamina and where, 
from (2-28), 

Ei 

*11 = 


4 - 4 , 


11 

1 - v no v 


12 21 


‘22 


<22 


<12 


' v 12 v 21 


Vi n E, 


12 22 


1 - v ] 2 v 21 


Q 66 ' G 12 

The five elastic constants which characterize the material are E^, E 2 2 » 

v 12* V 21 and G 12* ,our w ^ 11ch are independent. 

To impose isotropicity, we consider the limit of an infinite number 
of unidirectional bundles of fibers at equally spaced directions. Let j 
be the number of bundles having orientations between 0 and 0 + d0 and let 
m be the number of increments dO. Then mj = n and jtr * ndO so that 
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Carrying out the integration 



Since a bundle of fibers will have no transverse o» shear stiffness, 

Q 11 = E ll* Q 22 = Q 12 = Q 66 = 0 
so that 

E = — E 
L e 8 b ll 

Applying the rule of mixtures (ref. 13 shows this to be a good approximation) 
gives the stiffness of the composite as 

E C = I E 11 f V F X E M ^ “ V F^ 

where V F is the volume fraction of the filament. This expression was used 
in reference 14 and is used in this paper. 

There are other candidate expressions for and some of these will 
be briefly mentioned. I* the first, the limit of an infinite number of 
orthotropic lamina is considered instead of fiber bundles. In this case, 
transverse and shear properties cannot be neglected with the result that 


F = 1 A -1 

* T A 11 8 


E il . (3+2 v 12 ) E 22 , 1 c 

: + 8 T + 2 u 12 

1 - V] p V 21 * * 


Using the approximations (ref. 15) 


1 - v ]2 v 21 =1, v 12 = 1/4, G 12 = g E 22 


results in 


r = .3 r + 5. r 

e 8 t ll 8 to ? 


Applying the rule of mixtures to E^ 

E C = ! tE ll F V F + ^ + I 


72 
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where E 22 depends on Egg » ^ and Vp and may be determined from reference 
13 or 16. 

Other expressions for are obtained by considering the quasi -isotropic 
laminate subjected to axial loading. Then 


N x = A 11 e x + A 12 E y 


0 A 12 e x + A 22 e y 


or 


N x = (A 


12 x 

11 ‘ A 22 ; e x 


and we take the equivalent isotropic stiffness as 
E e = 7 (A 11 ' A 12* /A 22^ 

In the same manner as for A^, the following are obtained 

A 22 = A 11 = T ^ I + 1 ^12 + 2( W + I Q 22 ^ 


A 12 * E 8 ^11 + “ 2 Q 66 + 4 Q 12^ 

If bundles of fibers are considered and the transverse and shear stiffnesses 
are neglected, there results 



Applying the rule of mixtures gives 
E C = 3 E n p V F + ^ ^‘V 

If, on the other hand, layers of lamina are considered and the same 
approximations ar used as for E^ , there results 
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(g Ey) 4 g ^22^ ' F n + 


Q U 11 8 22 ' 


1 E + i E 

8 11 8 22 


mUa MO 
m ic re 


E u ■ E n. V F + ^ < 1 - V F> 


Calculations show that the various expressions for E c differ by about 10% 
for B/E and 3% for C/E. Since values of E,^ must *> e determined experimen- 
tally* it was decided not to use methods employing this quantity. 

For the composite density, the rule of mixtures is used 

P C = P F V F + P M ^ " V F^ 

It is assumed that composite strength is reduced by the same relative 
amount as the stiffness (ref. 11); thus 

F C = I F F V F + F M ^ " V F^ 

The properties of both the composite constituents and the composites 
themselves (with Vp = 0.5) are shown in Table Cl. Also shown are the 
properties of the aluminum alloy used in the study. Since the weight of 
buckling critical structure is approximately proportional to the quotient 
p^E, this value Is also given in the table, relative to aluminum. 
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APPENDIX D - LIFT AND DRAG ANALYSIS 
Lift Characteristics 

The lift characteristics, in terms of , were calculated using a thin- 

a 

wing paneling program developed by Ralph Carmichael at the Ames Research 
Center and based on the theory of reference 17. The analysis uses the 
method of aerodynamic influence coefficients, whereby the wing is divided 
into many small regions or panels with swept leading and trailing edges to 
conform to the actual wing geometry. Bilateral symmetry is assumed for all 
wings. The lifting pressure difference is assumed to be constant over each 
wing panel. The number of panels into which the wing is divided affects 
both the accuracy and the computer time required for the solution. As 
expected, the error decreases with increasing number of panels, although 
the convergence is rather slow. From a practical point of view, it is 
fortunate that answers with adequate accuracy for engineering analysis can 
be obtained with a relatively small number of panels, say 20-40. Since 
the computing time is roughly pror'rtional to the square of the number of 
panels, the use of a fine grid is justified only when a high degree of 
accuracy is required. In the analysis of cambered and twisted wings, however, 
it is usually necessary to use a large number of panels in order to adequately 
describe the wing camber surface. Since, in the use of this program for 
the analysis of advanced commercial transport configurations, it is currently 
not necessary to consider camber and twist, the wing lift characteristics 
can be adequately estimated with thirty-six panels. 
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Zero-Lift Drag 


The zero-lift drag (C n ), for each configuration was estimated by adding 

u 0 

the friction drag of the wing, fuselage, horizontal, and vertical tail, 
together with any applicable wave drag. The nacelle friction drag was included 
in the engine • 'crmance calculation. For configurations designed to 
cruise at speeds less than Mach 1, no wave drag was included. Also, no 
separation drag component was included in the zero-lift drag estimate. The 
estimated wing separation drag was included in the calculation of -nduced 
drag, as will be discussed subsequently. 

The friction drags for the aircraft components were computed using the 
following relations: 

The wing friction drag was estimated based on turbulent boundary layer, 
flat-plate skin friction and contains an empirical correction for thickness 
induced pressure fields. The basic equation is Frankl and Voishel's extension 
of Von Karman's mixing-length hypothesis to compressible flow, and the 
empirical correction for thickness-induced pressure fields is the result of a 
correlation in reference 18 of a large amount of data. The resulting equation 
is as follows: 



[i + 2 (i) ] ( !m) 

0.455 wi "9 

Oog )0 Re) 2 - 58 0 « 2 ) 


ufharo . 

1111% I V 9 


wing 

S WET 


= wing thickness ratio 

= wing wetted area, not including portion of wing 

2 

buried in the fuselage, ft 
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S 


REF 


Re 


Y 

M 


2 

= reference area, idealized wing planform area, ft 

= Reynold's Number based on the wing mean -aerodynamic- 
chord 

= ratio of specific heats for air 
= Mach number 


The fuselage friction drag was estimated with the same equation as used 
for the wing friction drag. However, the correction for thickness induced 
pressure fields is changed to the correct form for axisymmetric bodies: 


[ 1+2 (— ) . ] 
L v c'wmg J 


is replaced by 


+ 2 ^BODY-' 


where, (J-) B ody = ,fuse ^ a 9 e fineness ratio 

Also, for the fuselage friction drag the Reynold's Number used in the equation 
is based on body length and is the fuselage wetted area. 

The horizontal and vertical tail friction drags were estimated with the 
same equation as ised for the wing friction drag. In this case, the Reynold's 
Number is based on the appropriate mean-aerodynamic-chord and the thickness 
ration and wetted area used are those of the particular surface. 

No wave drag was included for the CVT configuration, which was designed 
to cruise at Mach 0.90. It was assumed that wave drag of the CVT with the 
supercritical wing would be negligible. Also, for the ATT configuration, 
which was designed to cruise at Mach 0.98, no wave drag was included. For 
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the ATT It was assumed that any wave drag would be negligible because of the 
supercritical wing and fuselage area ruling. 

The wave drag of the TST configuration, which was designed to cruise 
at Mach 1.15, was estimated as the theoretical wave drag of a axi syimietri c 
body with a Sears-Haack area distribution plus the wave drag of the 
horizontal and vertical tails, without any interference drag. 

The wave drag of a body with a Sears-Haack area distribution is the 
theoretical minimum for an axisymnetric body of a given length and volume. 
The form of the equation used for the Sears-Haack wave drag computation 
is as follows (ref. 19): 


C 



where , 


vol sh 

v°4h 



3 

volume of the Sears-Haack area distribution, ft 
length of the Sears-Haack area distribution, ft 
maximum diameter of the Sears-Haack area distribution, ft 


The equations used for estimating the wave drag of the horizontal and 
vertical tails are for values at Mach 1.0 and the Mach number for shock 


• ^ i. A A U*M 

Q U CttWi liiiun % 


+ U A ^ 

vUW* I til 
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Tho walno iiearl a+ Marh 1 u#ac nh+ainfiH 

• » » W t W • ...» ^ . . . W 


by linear interpolation between these two values. The equations for both 
fins are of the same form. Using the appropriate values of thickness ratio, 
exposed area, and sweep for that fin, the equations are as follows (ref. 20): 



- 4Q ~ 


At 


t 5 / 3 S fin ? 

M = 1.0, C n = 3.4 (£) * cos 2 

U WF c fin b ref 


A .. 

fin 


At M = secant A x . , 

fin* 

. 2 i S r . 

C D ' 6 (|) } 5^ 

U WF c fin e >ef 


where, (■£) = fin fineness ratio 

c fin 


S^n = fin ex P° se< ^ planform area, ft 


A 


= fin leading edge sweep, deg 


6 = compressibility factor ()M 2 - l|) 


1/2 


Induced Drag 

The induced drag of all configurations was estimated using an equation 
made up of the sum of three components: vortex drag, wave drag, and 
separation drag. The induced drag level is related to a parameter called 
K r def i ned by 


U D. 

2 

where C n = Induced drag coefficient and C, /.rAS, is the theoretical minimum 
u i L 

subsonic induced drag, with usual d^s efficiency factor equal io 1. 

The vortex drag was estimated to reduce Oswald's efficiency factor 

from 1.0 to 0.85, which corresponds to a K (the K factor for vortex 

drag) of 1.18. The lift Induced wave drag exponent was estimated at the 

theoretical minimum level given in references 21 and 22. This wave drag 

component was only inclided at speeds greater than Mach 1. The form of 
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this equation for conventional swept wings is as follows: 




ave 


where, 




w 7 ~~ j. x 3T2 
(tan a le + 1+x A ) 

= 1.0, val ue for wing shapes of nearly elliptic loading 
both spanwise and lengthwise, as defined in reference 22. 



wing leading edge sweep, deg 


A = wing taper ratio, tip chord/centerline chord 

p 

= wing aspect ratio, span/wing area 


The separation drag component was estimated from experimental data in 
reference 23. The separation drag estimate was included in an effort to 
realistically account for this drag component, which is difficult to compute 
theoretically. 

Figure D1 shows a comparison of the theoretical estima; of induced 
drag without including separation drag with some experimental data from 
reference 23. Based on this experimental data an empirical relation for 
estimating separation drag was derived. The resulting equation is as 
follows: 



C L } 3/2 

0 . 2 ) Cortex + *w ave 


1 . 0 ) 


where, = lift coefficient 




The expression for the induced drag coefficient resulting from combining 
the previously discussed components is as follows: 


'Dj ~ ^vortex + *wave + K sep^ 
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The induced drag of aircraft designed to cruise in the transonic speed 
regime is one of the most challenging aerodynamic quantities to predict. 
Realizing possible weakness in the estimation technique presented here, 
we have included several optional methods in the aerodynamics program. 

These optional estimation methods range from ideal subsonic induced drag 
with full leading edge suction (Oswald's efficiency factor equal to 1) 
to supersonic induced drag, corresponding to no leading edge suction 


(C D = C L tan a). The sensitivity to 
from these various methods is given in 


the induced drag estimate resulting 
the text of this report. 
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APPENDIX E - AREA RULING 

The area ruling program described herein was developed in response 
to the requirement for a simplified design tool to use as part of the 
TRANSYN-TST analysis of advanced technology transport aircraft. The purpose 
of the program is to calculate the fuselage cross-sectional area distribu- 
tion which, when combined with the wing and nacelle area distributions, 
will result in a given overall equivalent area distribution (a Sears- 
Haack area distribution in this study). The program is restricted to 
normal cross-sectional cuts and only includes the cross-sectional areas 
of the fuselage, wing, and engine nacelles. The fuselage area distribution 
is required for use in the fuselage structural weight calculation. Fuselage 
structural weight is dependent on the amount of area ruling required 
because the area ruling results in decreases in the fuselage diameter in 
the region where the fuselage bending moments are high. 

To maintain the same passengers capacity in an area-ruled fuselage as 
in an equivalent cylindrical fuselage, the floor area of each is made the 
same. This assumes that the efficiency of seat placement in an area-ruled 
configuration is the same as that for a cylindrical fuselage. Since this 
assumption is somewhat unrealistic, several area-ruled fuselages were 
examined to determine the actual number of seats that could be accommodated. 
Based on the results of this study, the floor area required for the area- 
ruled fuselage was increased slightly to allow for less efficient floor 
area utilization resulting from the varying fuselage diameter. For more 
detailed studies, it would be desirable to accurately compute the number 
of seats that could be accommodated for each fuselage area distribution. 

This would involve a station-by-station examination over the length of the 
fusel age . 
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In the TRANSYN-TST program, the area-ruling procedure is divided into 
several subroutines. These subroutines calculate the cross-sectional area 
distributions of the wing alone, of the engine nacelles (minus the flow- 
through area), and of the fuselage (allowing for the wing and nacelles). 

The area distribution of the tails is not included, but this area would be 
removed *rom the fuselage in region of low bending moment concentration, and 
this omission would have a very small effect on the calculated fuselage 
weight. 

For the calculation of the wing cross-sectional area distribution, 
the wing is assumed to be a constant airfoil section with a constant thickness 
ratio across the wing span. The chordwise thickness distribution is 
soecified bv the thickness ratin (t/c) at 12 stations over the chord of 

i » * 

the section. Between these stations the thickness is determined by straight 
line interpolation. The specification of the airfoil thickness distribution, 
combined with the geometric description of the wing planform (sweep and taper 
ratio), completely describes t’> thickness distribution of the wing. The 
wing cross-sectional area distribution is calculated at 11 equally spaced 
streamwise stations along the wing region of influence, from the leading 
edge of the wing centerline section to the trailing edge of the tip airfoil 
section. The cross-sectional area distributions for several wings of 
different sweeps and thickness ratios are shown in figure El. 

The geometry of the nacelles is determined by the propulsion system 
analysis (ref. 2) and assumes that the nacelle is axisymmetric. The nacelle 
exterior profile consists of a constant diameter engine section and circular 
arc radius calculated for tangency to the engine section while meeting the 
inlet and exit requirements. For the area distribution calculations, the 



interior shape of the flow-thru area is specified Dy a conical surface from 
the inlet lip to the nozzle lip. The effects of any exhaust gas expansion 
after leaving the nozzle are neglected in the are?.-ruling calculation. 

The total cross-sectional area distribution for multiple engines is 
dependent on the engine placement, and several subroutines have been 
developed for different types f engine installation* (4, and 8 wing- 
mounted engines, and 4 aft fuselage-mounted engines). 

When the engines are wing mounted the nacelle cross-section areas are 
computed for same streamwise stations as used in the wing area distribution 
calculation. The nacelle and wing cross-sectional area distributions are 
then combined directly when area ruling the fuselage. When the engines are 
mounted on the aft fuselage, the nacelle cross-sectional areas are computed 
at 11 equally spaced streamwise stations along the nacelle region of 
influence, and this nacelle area distribution is accounted for separately 
during the fuselage area ruling. 

To calculate the area distribution of the fuselage, the configuration 
geometry in the vicinity of the intersection of the wing and fuselage must 
be determined. For a low wing, the location of the wing was specified such 
that the bottom of a non-area ruled fu elage of the desired area distribu- 
tion coincided with the lowest point on the wing. To provide maximum 
headroom as the body was necked down, the centerline he fuselage was 
allowed to rise in the region above the wing. This displacement of the 
fuselage centerline was constrained so that the + op a- ;he area-ruKu 
fuselage was not higher than for the non-area ruled fuselage and the bottom 
of the fuselage remained in contact with the w’ng. The cross-sectional area 
of the combined wing and fuselage also includes the area of a fillet at the 
fuselage-wing junction. 
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The objecti .a of the fuselage area ruMng calculations was to determine 
the diameter distribution (d ARfi ) of the area-ruled body that results in a 
given total cross-sectional area distributio. In solving for this diameter 
distribution three limiting cases arise. For the case where the wing is 
completely included in the original body (shown in figure E2 as Case I), 
no fuselage modification is required and d ARfi is equal to d^g, the original 
fuselage diameter. When the wing cross-secticnal area is completely exterior 
to the fuselage (an aft cut, shown in figure E2 as Case II), the diameter 
of the area- ruled body d^g is calculated based on the original fuselage 
area minus the wing cross-sectional area. When the area-ruled fuselage is 
tangent to the wing surface the geometry is shown in figure E2 as Case III. 
All other foselage-wir ’’ntersections require a solution of the general case. 
Ihis is handled in the program by a numerical technique, iterating on o ARB 
until the required equivalent cross-sectional area is obtained. 

In order to maintain the required floor area for the area-ruled 
config ration, there is an additional iteration which adjusts the body length 
in ovder to obtain the required floor area. The currently used procedure 
for maintaining the correct fuselage capacity is shown schematically in 
figure E3. First, an input cylindrical body is sized on the basis of number 
of seats, seats abreast, seat pitch, seat width, number of aisles, and nose 
and afterbody fineness ratio (assuming Scars-Haack nose and afterbody area 
distributions). Second, the volume of this fuselage is calculated and a 
Srars-Haack body witn the same length and volume is computed. Third, the 
fineness ratio of \> Sears-Haack body is held constant while the fuselage 
is area ruled for the wing t *. nacelles and the body lennth is increased to 
maintain the same floor area. 




Table 1. Potential of Advanced Composites 
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Figure 2. Comparison of Actual and Estimated Empty Weights 
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Figure 3. Aerodynamic Estimation Methods 
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Figure 5. Aerodynamic Constraints 
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Figure 6. Nominal Configurations 
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Figure 7. Body Weight Characteristics 
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Figure 8a. Effect of Aspect Ratio on Weight Fractions 
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Figure 8b. Effect of Aspect Ratio on O.W.E., L/D max and Gross Weight 



RANGE = 2700 n.mi. 200 PASSENGERS 



Figure 9a. Effect of Wing Loading on Weight Fractions 
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Figure 9b. Effect of Wing Loading on O.W.E., L/D max and Gross Weight 
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Figure 10a. Effect of Thickness-to-Chord on Weight Fractions 
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Figure 10b. Effect of Thickness -to-Chord on O.W.E., L/D and Gross Weight 



RANGE = 2700 n.mi. 200 PASSENGERS 



CP 
< V 




Figure 11a. Effect of Sweep on Weight Fractions 
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Figure 11b. Effect of Sweep on O.W.E., L/D and Gross Weight 
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Figure 13. Effect of Aspect Ratio on O.W.E., L/D and Gross Weight 
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Figure 14. Effect of Wing Loading on O.W.E., L/D and Gross Weight 
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Figure 15. Effect of Thickness -to-Chord on O.W.E., L/D mav and Gross Weight 
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Figure 16. Effect of Sweep on O.W.E., L/D max and Gross Weight 
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Figure 17. heights of Optimized Configurations 
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Figure 18. Optimuir Aluminum Configurations 
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Figure 19. Optimum Carbon/Epoxy Configuration' 


RANGE = 2700 n.mi. 200 PASSENGERS 



S|0U>| ‘HOVOdddV^ 




XDU, * 1 0 


Figure 20. Effect of Wing Loading on Cj_ and Landing Speed 
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Figure 23. Effect of Landing Speed Constraint on Optimum Aluminum ATT Configuration 



RANGE = 2700 n. mi. 200 PASSENGERS 



CM :9 

CO 

5= O 

^ O 


■*> c 

— • JO o 00 


-~sm 

— CD 

£ 

if) 

oo ro 

— 


CVI_ £ (J) 

CO \ 8 O CO g 

oi £ O <g CM ^ 

Q fO If) GO 

O CVJ CD 

— ro — 



I- 

o 

LU 

CL 

CO 

< 



O 




O 
h- 

O UJ 

$ CO 

5 s => 
^ cr 

. o 

CO UJ ^ 

. o 
£ 


o 

LU 

£ 


CO 

o 

a: _ 

o o d 



Figure 24. Effect of Landing Speed Constraint on Optimum Carbon/Epoxy JST Configuration 
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Figure B2. Wing Coordinate System and Structural Concept 
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Figure E2. Fuselage/Wing Intersection Geometry Limiting Cases 




Figure E3. Area Ruling Procedure 




